function  [qp,qw,b1,b2,b3]=quadrule
p=[0.33333333333333   0.33333333333333   0.03290633138892
   0.02078002585399   0.48960998707301   0.01033073189127
   0.48960998707301   0.48960998707301   0.01033073189127
   0.48960998707301   0.02078002585399   0.01033073189127
   0.09092621460422   0.45453689269789   0.02238724726302
   0.45453689269789   0.45453689269789   0.02238724726302
   0.45453689269789   0.09092621460422   0.02238724726302
   0.19716663870114   0.40141668064943   0.03026612586947
   0.40141668064943   0.40141668064943   0.03026612586947
   0.40141668064943   0.19716663870114   0.03026612586947
   0.48889669119381   0.25555165440310   0.03049096780220
   0.25555165440310   0.25555165440310   0.03049096780220
   0.25555165440310   0.48889669119381   0.03049096780220
   0.64584411569574   0.17707794215213   0.02415921274164
   0.17707794215213   0.17707794215213   0.02415921274164
   0.17707794215213   0.64584411569574   0.02415921274164
   0.77987789354410   0.11006105322795   0.01605080358680
   0.11006105322795   0.11006105322795   0.01605080358680
   0.11006105322795   0.77987789354410   0.01605080358680
   0.88894275149632   0.05552862425184   0.00808458026178
   0.05552862425184   0.05552862425184   0.00808458026178
   0.05552862425184   0.88894275149632   0.00808458026178
   0.97475627244554   0.01262186377723   0.00207936202748
   0.01262186377723   0.01262186377723   0.00207936202748
   0.01262186377723   0.97475627244554   0.00207936202748
   0.00361141784841   0.39575478735694   0.00388487690498
   0.39575478735694   0.60063379479464   0.00388487690498
   0.60063379479464   0.00361141784841   0.00388487690498
   0.39575478735694   0.00361141784841   0.00388487690498
   0.60063379479464   0.39575478735694   0.00388487690498
   0.00361141784841   0.60063379479464   0.00388487690498
   0.13446675453078   0.30792998388044   0.02557416061202
   0.30792998388044   0.55760326158878   0.02557416061202
   0.55760326158878   0.13446675453078   0.02557416061202
   0.30792998388044   0.13446675453078   0.02557416061202
   0.55760326158878   0.30792998388044   0.02557416061202
   0.13446675453078   0.55760326158878   0.02557416061202
   0.01444602577612   0.26456694840652   0.00888090357334
   0.26456694840652   0.72098702581736   0.00888090357334
   0.72098702581736   0.01444602577612   0.00888090357334
   0.26456694840652   0.01444602577612   0.00888090357334
   0.72098702581736   0.26456694840652   0.00888090357334
   0.01444602577612   0.72098702581736   0.00888090357334
   0.04693357883818   0.35853935220595   0.01612454676173
   0.35853935220595   0.59452706895587   0.01612454676173
   0.59452706895587   0.04693357883818   0.01612454676173
   0.35853935220595   0.04693357883818   0.01612454676173
   0.59452706895587   0.35853935220595   0.01612454676173
   0.04693357883818   0.59452706895587   0.01612454676173
   0.00286112035057   0.15780740596859   0.00249194181749
   0.15780740596859   0.83933147368084   0.00249194181749
   0.83933147368084   0.00286112035057   0.00249194181749
   0.15780740596859   0.00286112035057   0.00249194181749
   0.83933147368084   0.15780740596859   0.00249194181749
   0.00286112035057   0.83933147368084   0.00249194181749
   0.22386142409792   0.07505059697591   0.01824284011895
   0.07505059697591   0.70108797892617   0.01824284011895
   0.70108797892617   0.22386142409792   0.01824284011895
   0.07505059697591   0.22386142409792   0.01824284011895
   0.70108797892617   0.07505059697591   0.01824284011895
   0.22386142409792   0.70108797892617   0.01824284011895
   0.03464707481676   0.14242160111338   0.01025856373620
   0.14242160111338   0.82293132406986   0.01025856373620
   0.82293132406986   0.03464707481676   0.01025856373620
   0.14242160111338   0.03464707481676   0.01025856373620
   0.82293132406986   0.14242160111338   0.01025856373620
   0.03464707481676   0.82293132406986   0.01025856373620
   0.01016111929628   0.06549462808294   0.00379992885530
   0.06549462808294   0.92434425262078   0.00379992885530
   0.92434425262078   0.01016111929628   0.00379992885530
   0.06549462808294   0.01016111929628   0.00379992885530
   0.92434425262078   0.06549462808294   0.00379992885530
   0.01016111929628   0.92434425262078   0.00379992885530];
qp=p(:,1:2);
qw=p(:,3);
[b1,b2,b3]=barycenter(p(:,1),p(:,2));
% x0=[0 1 0];
% y0=[0 0 1];
% W=diag(w);
% [b1 b2 b3]=barycenter(x1,y1,x0,y0);
% mat=vdm23(d,b1,b2,b3);
% G=mat'*W*mat
% max(max(abs(build(d)-G)))
